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Abstract 

We have used the event generator LUCIAE to analyse NA35 data of 
the p and A yields, the ratio A/p, and the transverse mass distributions of p 
and A in pp and central sulphur- nucleus collisions at 200 A GeV. The NA35 
data could be reproduced reasonably if one assumes that the s quark pair 
suppression factor and concerns in nucleus-nucleus collisions are larger 
than the nucleon-nucleon collision. It seems to indicate that NA35 data 
might imply the reduction of strangeness suppression in ultrarelativistic 
nucleus-nucleus collisions comparing to the nucleon-nucleon collision at 
the same energy. However, the ratio A/p approaching to unity in AA 
collisions comparing to the pp collision does not necessarily mean a flavor 
symmetry since hadronic rescattering plays a role as well. 
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Strangeness production is expected to be a powerful probe for the mech- 
anism of nucleus- nucleus collisions. Strangeness enhancement, the increased 
strangeness particle production in nucleus-nucleus collisions comparing to the 
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nucleon-nucleon collision, is predicted to be a sensitive signature of the QGP 
formation in the ultrarelativistic nucleus-nucleus collisions [1]. 

The first experimental results of the enhanced production of strange parti- 
cles in nucleus-nucleus collisions at 200A GcV incident energy was reported six 
years ago [2]. Later on this enhancement was confirmed by more and more 
experiments [3-7]. 

Strange antibaryon production, in particular, might bring messages of the 
equilibrium and flavour symmetry of quarks in ultrarelativistic nucleus-nucleus 
collisions. It has been estimated [8] that if there is QGP formed in the ultra- 
relativistic nucleus- nucleus collisions, strange (antistrange) quark pair might be 
copiously reproduced, resulting in an approximate flavour symmetry among u, 
d and s quarks [9]. Thus the ratio of A/p should approach unity, since A is 
composed of uds and p is udd. 

Recently published NA35 data of antibaryon production in sulphur-nucleus 
collisions at 200A GeV [7] did really observe that the ratios A /p are all ap- 
proaching, even over, unity and far exceeding the corresponding value of ~0.25 
in the nucleon-nucleon collision at the same energy. Although the NA35 data 
of the A/p has been fairly reproduced by RQMD 1.08 [10,7] with fusion of 
overlpping strings into a color rope and with hadronic rescattering, the data 
of transverse mass distributions of p and A have not been explained yet, the 
physics behind the data have not been exposed especially. 

Based on the idea that the strangeness enhancement in ultrarelativistic nucleus- 
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nucleus collision compared to the nucleon-nucleon collision at the same energy 
could be investigated via the reduction of s quark suppression [11-14], in this 
letter the event generator LUCIAE [15] is used to analyse the NA35 data and 
to explore the physics behind the data. 

LUCIAE was updated based on FRITIOF 7.02 [16] by taking into account the 
collective string interaction [17-18] during the emission of gluon bremsstrahlung 
and the rescattering of produced particles [19]. One knows that in FRITIOF 
7.02 [16], two colliding hadrons are excited due to the longitudinal momentum 
transfer and/or a Rutherford parton scattering. The highly excited states will 
emit bremsstrahlung gluons according to the soft radiation model. The deex- 
cited states are then treated as Lund strings allowing to decay into the final 
hadronic state due to the Lund fragmentation scheme. However, in the ultra- 
relativistic nucleus-nucleus collisions there are generally many excited strings 
formed close by each other. These strings would behave like vortex lines in a 
color superconducting QCD vacuum and will interact with each other (as the 
repulsive interaction suffered by the "ordinary" vortex lines in a type II super- 
conductor). This kind of collective interaction among strings, not included in 
FRITIOF 7.02, is depicted by the firecracker model [15,18] dealing with the large 
p t gluon (firecracker gluon) production from the collective interaction among 
strings. The firecracker gluon will work as a gluon-kink excitation on a string, 
which tends to shorten the longitudinal size of the string and then brings about 
the increasing of the string tension (thus it should be regarded as an effective 
string tension). In firecracker model it is assumed that the groups of neighbor- 
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ing strings might form interacting quantum states, the large common energy 
density (corresponding to the collective interaction among strings) might then 
affect the emission of gluonic bremsstrahlung [17-18]. 

A rescattering model has been developed to describe the reinteraction of pro- 
duced particles, from FRITIOF event generator, with each other and with the 
participant and the spectator nucleons [19]. In the model, the produced parti- 
cles and the participant (wounded) nucleons, from FRITIOF event generator, 
are distributed randomly in the geometrical overlapping region between the 
projectile and the target nuclei under a given impact parameter. The target 
(projectile) spectator nucleons are distributed randomly outside the overlap- 
ping region and inside the target (projectile) sphere. A rescattering cascade 
process has evolved since then, cf. Ref. [15,19] for the detail. The considered 
inelastic reactions, concerning strangeness, are here cataloged into: 



7T7T # kk; 




nN # kY, 


ttN # kY; 


ttY # kE, 


ttY # kE; 


kN^irY , 


kN ^ nY; 


kY # irE, 


kY # ttE; 


kN ^ jfcS, 


kN ^ kE; 


nE # k£l~, 


nE # kn~: 


kE # irQ~, 


kE ^ irQ,~: 


NN annihilation; 
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YN annihilation; 



where Y refers to the A or E and S refers to the S~ or S°. There are 299 
inelastic reactions involved altogether. As the reactions introduced above do not 
make up the full inelastic cross section, the remainder is again treated as elastic 
scattering [19]. The cross section of -kit — > kk is taken to be 2.0 mb as usual 
[10]. The isospin averaged parametrization of Ref. [9] is adopted for the cross 
sections of the reactions ttN — > kY and for the other strange quark production 
reactions. Of course, the difference in threshold energy among reactions is taken 
into account. Following Ref. [9], the cross section of strange quark exchange 
reaction, kN — > nY for instance, is assumed to be equal to ten times the value 
of the cross section of the strangeness production reaction. As for the cross 
section of the reverse reaction, the detailed balance assumption [20] is required. 

The cross sections of the inelastic reactions given by the isospin averaged 
parameterization formulas of [9] for the irN — > kY decrease exponentially with 
the CMS energy of the two colliding particles. But we know that the total 
inelastic cross section of the ttN is approximately energy -independent, which 
means that more inelastic channels would occur like n~ + p — > K*° + A — > K + 
7T + A and n~ + p -> K*° + S° -> K + tt + S° etc., when the CMS energy of the 
two colliding particles is increasing. In some sense the reactions we listed above 
should be looked upon as 'representative channels' of certain types of reactions. 
Therefore instead of applying energy-dependent inelastic cross section, we could 
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alternatively give a constant 'effective cross section' to the reaction irN — > kY 
[15,19]. 

In JETSET routine, which runs together with LUCIAE event generator, there 
are model parameters parj(2) (or 's') and parj(3) which are responsible for the 
s quark suppression and related to the effective string tension. V refers to the 
suppression of s quark pair production in the color field compared to u or d pair 
production. parj(3) is the extra suppression of strange diquark production com- 
pared to the normal suppression of strange quark pair. Besides V and parj(3) 
there is parj(l), which stands for the suppression of diquark- antidiquark pair 
production in the color field compared to the quark-antiquark pair production 
and is related to the effective string tension as well. Originally, in the LUND 
fragmentation scheme 's' was assumed to be a 'constant' and this assumption 
was confirmed by the e + e~ physics from the low energies to Z° energy. However, 
there are experimental facts that this parameter is energy dependent when the 
fragmentation scheme is applied to the Deep Inelastic Scattering (DIS) and hh 
collisions. In DIS experiments, 's' is seen to be rising up with the increase of 
energy from 0.15 at V5=5 GeV to 0.35 at ^=20 GeV [21]. An energy depen- 
dent 's' for hh collisions has been known for many years. It varies from 0.2 at 
the ISR energy to about 0.4 at v / 5= 1-8 TeV [21]. In addition, this parameter 
is also observed strongly phase space dependence in DIS and hh collisions. The 
corresponding value of 's' runs from 0.15 to 0.55 with a mean value close to 0.3 
[22] . Thus the idea of the reduction of s quark suppression in nucleus-nucleus 
collisions comparing to the nucleon-nucleon collision is possible to be executed 
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via changing the 's' and concerning parameters in JETSET routine. 

From LUND string model point of view s quark suppression factor could be 
related to the effective string tension. It is reasonable to expected that the 
more violent collision, the stronger collective interaction among strings, the 
larger effective string tension and then the larger s quark suppression factor. 
Generally speaking, a nucleus-nucleus collision is more violent than the nucleon- 
nucleon collision at the same interaction energy. That might be the reason of 
the reduction of s quark suppression in nucleus- nucleus collisions comparing to 
the nucleon-nucleon collision at the same energy. 

The purpose of this letter is to explore the physics behind the NA35 data 
and not to fit the data as good as possible. We fix a set of somewhat larger 
parameters (relative to the defaults: 's'=0.3, parj(3)=0.4 and parj(l)=0.1) of 
's'=0.4, parj(3)=0.5333 and parj(l)=0.1333 (referred to as parameter set 1, 
later on) in the calculations of nucleus- nucleus collisions and a set of somewhat 
smaller parameters of 's'=0.2, parj(3)=0.2666 and parj(l)=0. 06666 ( referred to 
as parameter set 2) in the calculation of the nucleon-nucleon collision . These 
results are given in table 1 and figure 1 and obtained from the average over 
10 5 generated events for pp, 2000 events for A and p yield in AA, and 3000 
for m t distribution in AA. The y and p t acceptances are set to be the same 
as the experiment [7], correspondingly. Of course, the two sets of parameters 
above imply that one makes an assumption there that those three parameters 
are linearly proportional to the effective string tension in their responsibilities 
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to the final hadronic production. 

Tab. 1 gives the rapidity densities of p, A, and h~ and the ratio ~h.jp in p 
+ p and S + S, S + Ag, and S + Au central collisions at 200 A GeV. The 
rapidity acceptance is 3 < y < 4. The results of LUCIAE are comparable 
with the corresponding NA35 data. The LUCIAE results of A jp seem going 
down monotonously from S + S to S + Au, which might attribute to the fact 
that one did not consider the possible difference of the reduction of strangeness 
suppression among them. The results of LUCIAE and RQMD are comparable 
with each other, since both of them do have contained the collective interaction 
among strings and the hadronic rescattering [23]. The hadronic rescattering, in 
LUCIAE and RQMD, is expected to play a similar role in the final state dis- 
tributions though details of the hadronic rescattering are not the same in these 
two generators. In RQMD the effect of highly dense strings in ultrarelativistic 
nucleus-nucleus collisions is considered via the fusion of strings into a color rope, 
which brings the enhanced production of strange antibaryon. In LUCIAE the 
highly dense strings are considered by the firecracker model [15,18], where the 
effect of highly dense strings is depicted analytically only in the emission of gluon 
bremsstrahlung. It leaves the possibility to enlarge the s quark pair suppression 
factor and concerns in JETSET routine to bring the enhanced production of 
strange antibaryon. 

Fig. 1 gives p and A transverse mass distributions in central S + S (upper 
frame), S + Ag (middle frame), and S + Au (lower frame) reactions at 200 A 
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GcV, respectively. In this figure, the open squares and open circles are the 
NA35 data of p and A, respectively and the corresponding results of LUCIAE 
are given by full squares and full circles. The solid lines are the exponential fits 
of the NA35 data, cf. Ref. 7 for details. One sees from this figure that the 
agreement between the NA35 data and the results of LUCIAE is fair for p and 
reasonably good for A. 

Table 2 gives the average yield (in full phase space) of A, A, K + and K°, in p+p 
and central S + S and S + Ag collisions at 200A GeV. In this table the data were 
taken from [5] besides p+p, which was taken from [24] . One sees from this table 
that although we are here aiming at analysing strange antibaryon production 
data, the strange particle production data are also reproduced reasonably good 
at the same time. 

In order to distinguish the role of the s quark pair suppression factor and 
concerns from the hadronic rescattering and the role of V and parj(3) from 
parj(l) one calculates the results of table 3 and figure 2. In table 3 one compares 
the results of the rapidity densities of p, A, and h~ and the ratio A/p in S + S 
central collisions at 200A GeV calculated by using different parameters and with 
or without rescattering: 'LUCIAE 1': parameter set 2 and without rescattering; 
'LUCIAE 2': parameter set 2 and with rescattering; 'LUCIAE 3': parameter 
set 1 but the value of parj(l) is changed to 0.06666 and with rescattering; and 
'LUCIAE 4': parameter set 1 and with rescattering. One sees from this table 
that the 'LUCIAE 1' result of A/p is close to the corresponding result in the 
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pp collision, as it should be. For results of 'LUCIAE 2', A yield is close to 
the corresponding result of 'LUCIAE 1' but p yield is lower than 'LUCIAE 
1', which brings about the larger ratio A/p in 'LUCIAE 2' than in 'LUCIAE 
1'. Here one knows that the hadronic rescattering seems to play nearly null 
role for A yield, since A production in rescattering is mainly via irN ^ kY 
and kN ^ irY 7 that is nearly canceled by the corresponding inverse reactions 
and YN annihilations especially. On the contrary, the rescattering tends to 
reduce p multiplicity through the pp annihilation. However, relying on hadronic 
rescattering only is not possible to have the ratio A/p approaching to unity. 
Furthermore, comparing the results of 'LUCIAE 3' with 'LUCIAE 4' one knows 
that although by increasing parj(l) the yield of p and A are increased the ratio 
A/p is hardly affected. 

Figure 2 gives the transverse mass distributions of A (upper frame) and p 
(lower frame) in S + S reaction at 200 A GeV, respectively. The full triangles, 
the open squares, and the full circles are calculated individually for the case 
1: parameter set 2 and with rescattering; the case 2: parameter set 1 and 
without rescattering; and the case 3: parameter set 1 and with rescattering. 
The rapidity acceptance are all set to be 1 < y < 3. One knows from this 
figure again that the hadronic rescattering seems to play nearly null role for A 
production, but hadronic rescattering tends to reduce p multiplicity, at lower 
m t region especially, through the pp annihilation. 

In summary, we have roughly reproduced the NA35 data of p, A yields , the 
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corresponding ratio A/p, and the transverse mass distributions of p and A in pp 
and central S + S, S + Ag, and S + Au collisions at 200 A GcV, using event gen- 
erator LUCIAE and via increasing s quark pair suppression factor and concerns 
in AA collisions. It seems to be true that the NA35 data [5,7] imply the re- 
duction of strangeness suppression in ultrarelativistic nucleus-nucleus collisions 
comparing to the pp collision at the same energy. Although the ratio A/p ap- 
proaching to unity is the result of reduction of s quark suppression, it does not 
necessarily mean a flavor symmetry, since reproducing the NA35 data does not 
require the 's' value should be equal to one and the hadronic rescattering also 
plays a role in enlarging the ratio A/p. However, it is absolutely needed to have 
a further study for the microscopic mechanism of the reduction of strangeness 
suppression before making a conclusion using strangeness enhancement as a 
signal of QGP. 
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Figure Captions 



Fig. 1 Transverse mass distributions of p (3 < y < 4) and A 
produced in central S+S (upper frame, 1 < y < 3), S+Ag (middle 
frame, 1 < y < 3), and S+Au (lower frame, 3 < y < 5) collisions at 
200A GeV. The open squares and open circles are the NA35 data 
of p and A, respectively and the full squares and full circles are the 
corresponding results of LUCIAE. 

Fig. 2 Transverse mass distributions (1 < y < 3) of A (upper 
frame) and p (lower frame) in S+S reaction at 200A GeV. The full 
triangles, the open squares, and the full circles are calculated indi- 
vidually for the case 1: parameter set 2 fs'=0.2, parj(3)=0.2666 and 
parj(l)=0. 06666) and with rescattering; the case 2: parameter set 1 
('s'— 0.4, parj(3)=0.5333, parj(l)=0.1333) and without rescattering; 
and the case 3: parameter set 1 and with rescattering. 
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Table 1 . Rapidity densities (3 < y < 4) of p, A and h 



produccc 


1 in p+p anc 


1 central sulphur-nucleus col 


lisions at 


200A GeV 


reaction 




V 
F 


A 


A/p 


h~ 


P + P 


data 


0.02±0.02 


005±0 002 


0.25±0.1 


0.74±0.04 




LUCIAE 


0.017 


0.0035 


0.21 


0.60 




RQMD 


0.015 


0.005 


0.3 




S + S 


data 


0.4±0.1 


0.76±0.16 




25±1 


central 


LUCIAE 


0.65 


0.66 


1.03 


23.7 




RQMD 


0.7 


0.75 


1.1 




S + Ag 


data 


0.6±0.2 


0.75±0.19 


1 3 +0 - 7 


40±2 


central 


LUCIAE 


1.08 


0.99 


0.91 


39.3 




RQMD 


1.0 


0.9 


0.9 




S + Au 


data 


0.7±0.2 


0.75±0.1 


1 1+°- 4 


47±5 


central 


LUCIAE 


1.09 


0.91 


0.84 


42.7 




RQMD 


1.4 


1.2 


0.9 




Table 2. Average yield (in full phase space) of A, A, K + and K° 


produced 


in p+p and central sulph 


ur-nucleus collisions at 200A GeV 


reaction 




A 


A 


K+ 




P + P 


data 


0.096±0.015 


0.013±0.005 




0.17±0.01 




LUCIAE 


0.10 


0.011 


0.22 


0.16 


S + S 


data 


9.4±1.0 


2.2±0.4 




10.5±1.7 


central 


LUCIAE 


8.0 


1.9 


11.9 


9.83 


S + Ag 


data 


15.2±1.2 


2.6±0.3 




15.5±1.5 


central 


LUCIAE 


13.9 


3.1 


21.3 


17.6 
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Table 3. Rapidity densities (3 < y < 4) of p, A and h 



produced in p+p and center S+S collision at 200A GeV 


reaction 




P 


A 


A/p 


h" 


P + P 


data 


0.02±0.02 


0.005±0.002 


0.25±0.1 0.74±0.04 




LUCIAE 


0.017 


0.0035 


0.21 


0.60 




RQMD 


0.015 


0.005 


0.3 




S + S 


data 


0.4±0.1 


0.76±0.16 


L - y -o.e 


25±1 


central 


LUCIAE l a 


0.62 


0.17 


0.28 


25.2 




LUCIAE 2 b 


0.42 


0.17 


0.40 


24.8 




LUCIAE 3 C 


0.36 


0.38 


1.05 


24.7 




LUCIAE 4 d 


0.65 


0.66 


1.03 


23.7 



a. 's'=0.2, parj(3)=0.2666, parj(l)=0. 06666 without rcscattcring 

b. 's'=0.2, parj(3)=0.2666, parj(l)=0.06666 with rescattering 

c. 's'=0.4, parj(3)=0.5333, parj(l)=0.06666 with rescattering 

d. 's'=0.4, parj(3)=0.5333, parj(l)=0.1333 with rescattering 
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